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Abstract. We investigate the properties of the zeros of the eigenfunctions on quan- 
tum graphs (metric graphs with a Schrodinger-type differential operator) . Using tools 
such as scattering approach and eigenvalue interlacing inequalities we derive several 
formulas relating the number of the zeros of the n-th eigenfunction to the spectrum 
of the graph and of some of its subgraphs. In a special case of the so-called dihedral 
graph we prove an explicit formula that only uses the lengths of the edges, entirely 
bypassing the information about the graph's eigenvalues. The results are explained 
from the point of view of the dynamics of zeros of the solutions to the scattering 
problem. 



1. Introduction 

Spectral properties of differential operators on graphs have recently arisen as models 
for such diverse areas of research as quantum chaos, photonic crystals, quantum wires 
and nanostructures. We refer the interested reader to the reviews [H [2] as well as to 
collections of recent results [31 H]. As a part of this research program, the study of 
eigenfunctions, and in particular, their nodal domains is an exciting and rapidly de- 
veloping research direction. It is an extension to graphs of the investigations of nodal 
domains on manifolds, which started already in the 19th century by the pioneering work 
of Chladni on the nodal structures of vibrating plates. Counting nodal domains started 
with Sturm's oscillation theorem which states that a vibrating string is divided into 
exactly n nodal intervals by the zeros of its n-th vibrational mode. In an attempt to 
generalize Sturm's theorem to manifolds in more than one dimension, Courant formu- 
lated his nodal domains theorem for vibrating membranes, which bounds the number of 
nodal domains of the n-th. eigenfunction by n [5]. Pleijel has shown later that Courant's 
bound can be realized only for finitely many eigenfunctions [5J. The study of nodal do- 
mains counts was revived after Blum et al have shown that nodal count statistics can 
be used as a criterion for quantum chaos [7J. A subsequent paper by Bogomolny and 
Schmit illuminated the fascinating connection between nodal statistics and percolation 
theory [8]. A recent paper by Nazarov and Sodin addresses the counting of nodal do- 
mains of eigenfunctions of the Laplacian on S> 2 [9]. They prove that on average, the 
number of nodal domains increases linearly with n, and the variance about the mean 
is bounded. At the same time, it was shown that the nodal sequence - the sequence of 
numbers of nodal domains ordered by the corresponding spectral parameters - stores 
geometrical information about the domain [10] . Moreover, there is a growing body of 
numerical and theoretical evidence which shows that the nodal sequence can be used 
to distinguish between isospectral manifolds [TT | [T2 | H3]. 

As far as counting nodal domains on graphs is concerned, it has been shown that trees 
behave as one-dimensional manifolds, and the analogue of Sturm's oscillation theory 
applies [H I [T5 | [T6 | ITT] . as long as the eigenfunction does not vanish at any vertex. Thus, 
denoting by v n the number of nodal domains of the n'th eigenfunction, one has u n — n 
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for tree graphs. Courant's theorem applies for the eigenf unctions on a generic graph: 
v n < n i [18J- It should be noted that there is a correction due to multiplicity of the n-th 
eigenvalue and the upper bound becomes n+m—1, where m is the multiplicity [19J. In 
addition, a lower bound for the number of nodal domains was discovered recently. It is 
shown in [2D] that the nodal domains count of the n-th eigenfunction has no less than 
n — (3 nodal domains, where is the number of independent cycles in the graph. Again, 
this result is valid for generic eigenfunctions, namely, the eigenfunction has no zero 
entries on the vertices and belongs to a simple eigenvalue. In a few cases, the nodal 
counts of isospectral quantum graphs were shown to be different, and thus provided 
further support to the conjecture that nodal count resolves isospectrality [21]. A recent 
review entitled "Nodal domains on graphs - How to count them and why?" [22] provides 
a detailed answer to the question which appears in its title (as it was known when the 
article was written). In particular, this manuscript contains a numerically established 
formula for the nodal count of a specific quantum graph, expressed in terms of the 
lengths of its edges. This was the first, and to this date the only, explicit nodal count 
formula for a non-trivial graph and in this manuscript we succeed in rigorously proving 
it. 

This leads us to focus on the study of nodal domains on quantum graphs from 
a new point of view. Namely, we shall show that one can count the number of nodal 
domains by using scattering data obtained by attaching semi-infinite leads to the graph. 
Scattering on graphs was proposed as a paradigm for chaotic scattering in [25| [21] with 
new applications and further developments in the field described in [2o~| |2U| [27]. The 
work presented here is based on the concepts and ideas developed in these studies. 

The paper is organized in the following way. The current section provides the nec- 
essary definitions and background from the theory of quantum graphs. The conversion 
of finite graphs to scattering systems by adding leads will be discussed in the next sec- 
tion and the expression for the scattering matrix will be derived and studied in detail. 
The connection of the scattering data with nodal domains and the counting methods 
it yields will be presented in section [3J Section H] applies the above counting methods 
in order to derive a formula for the nodal count of graphs with disjoint cycles. This 
formula relates the nodal count to the spectra of the graph and some of its subgraphs. 
Thus, information about the eigenfunctions is exclusively obtained from the eigenvalue 
spectrum. The last section relates the different ways of counting and discusses possible 
future developments. 

1.1. Quantum graphs. In this section we describe the quantum graph which is a 
metric graph with a Shrodinger-type self-adjoint operator defined on it. Let T = (V,£) 
be a connected graph with vertices V = {vj} and edges £ = {e,}. The sets V and 8 
are required to be finite. 

We are interested in metric graphs, i.e. the edges of V are 1-dimensional segments 
with finite positive lengths {L e } eg £. On the edge e = (u, v) we use two coordinates, x CjV 
and x eu . The coordinate x ev measures the distance along the edge starting from the 
vertex v ; x CtU is defined similarly. The two coordinates are connected by x e ^ v + x e ^ u = L e . 
Sometimes, when the precise nature of the coordinate is unimportant, we will simply 
write x e or even x. 

A metric graph becomes quantum after being equipped with an additional structure: 
assignment of a self-adjoint differential operator. This operator will be often called the 
Hamiltonian. In this paper we study the zeros of the eigenfunctions of the negative 
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where V(x) is a potential. Note that the value of a function or the second derivative 
of a function at a point on the edge is well-defined, thus it is not important which 
coordinate, x StV or x e n is used. This is in contrast to the first derivative which changes 
sign according to the direction of the chosen coordinate. 

To complete the definition of the operator we need to specify its domain. 

Definition 1.1. We denote by H 2 (T) the space 



which consists of the functions / on T that on each edge e belong to the Sobolev space 
H 2 (e). The restriction of / to the edge e is denoted by f e . The norm in the space 



Note that in the definition of H 2 (T) the smoothness is enforced along edges only, 
without any junction conditions at the vertices at all. However, the standard Sobolev 
trace theorem (e.g., [28]) implies that each function f e and its first derivative have well- 
defined values at the endpoints of the edge e. Since the direction is important for the 
first derivative, we will henceforth adopt the convention that, at an end- vertex of an 
edge e, the derivative is calculate into the edge and away from the vertex. That is the 
coordinate x is chosen so that the vertex corresponds to x — 0. 

To complete the definition of the operator we need to specify its domain. All condi- 
tions that lead to the operator (11. ip being self-adjoint have been classified in [29ll30tl31]. 
We will only be interested in the so-called extended 5-type conditions, since they are 
the only conditions that guarantee continuity of the eigenfunctions, something that is 
essential if one wants to study changes of sign of the said eigenfunctions. 

Definition 1.2. The domain % of the operator (II .2|) consists of the functions / G H 2 {T) 
such that 

(1) / is continuous on every vertex: 



# 2 (r) :=©# 2 (e), 



H 2 



(r) is 






(2) 



for every vertex v and edges e\ and e2 that have v as an endpoint. 
the derivatives of / at each vertex v satisfy 



(1.4) 




where S v is the set of edges incident to v. 
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Sometimes the condition (jl.4p is written in a more robust form 

(1.5) cos(7„) Y J^T^ = sin Mf(v), 

which is also meaningful for infinite values of a v = tan(7„). Henceforth we will under- 
stand a v = oo as the Dirichlet condition f(v) =0. The case a v = is often referred to 
as the Neumann-Kirchhoff condition. 

Finally, we will assume that the potential V(x) is bounded and piecewise continuous. 
To summarize our discussion, the operator (II. 2p with the domain T-L is self-adjoint for 
any choice of real a v . Since we only consider compact graphs, the spectrum is real, 
discrete and with no accumulation points. We will slightly abuse notation and denote 
by cr(r) the spectrum of an operator H defined on the graph T. It will be clear from 
the context which operator H we mean and what are the vertex conditions. 

The eigenvalues A G er(r) satisfy the equation 

(1.6) --^ + V(x)f(x) = \f(x). 

It can be shown that under the conditions specified above the operator H is bounded 
from below [31j. Thus we can number the eigenvalues in the ascending order, starting 
with 1. Sometimes we abuse the notation further and also call k, such that A = k 2 , 
an eigenvalue of the graph T. This also should lead to no confusion since, with the 
conditions Re (A;) > 0, lm(k) > 0, the relation between k and A is bijective. 

1.2. Nodal count. The main purpose of this article is to investigate the number of 
zeros and the number of nodal domains of the eigenfunctions of a connected quantum 
graph. We aim to give formulas linking these quantities to the geometry of the graphs 
and to the eigenvalues of the graph and its subgraphs, but avoiding any reference to 
the values of the eigenfunctions themselves. 

The number of internal zeros or nodal points of the function / will be denoted by n{f). 
We will use the shorthand /i n to denote /i(/ n ) where /„ is the n-th eigenfunction of the 
graph in question. The sequence {/x n } will be called the nodal point count sequence. A 
positive (negative) domain with respect to / is a maximal connected subset in V where 
/ is positive (correspondingly, negative). The total number of positive and negative 
domains will be called the nodal domain count of / and denoted by v(f). Similarly to 
fi n , we use v n as a short-hand for u(f n ) and refer to {u n } as the nodal domain count 
sequence. 

The two quantities fi n and v n are closely related, although, due to the graph topology, 
the relationship is more complex than on a line, where v = \i + 1. Namely, one can 
easily establish the bound 

(1.7) H-Pr + I<v<n + 1, 

where flr is the cyclomatic number of T. The cyclomatic number can be computed as 

(1.8) Pr=\£\-\V\ + 1. 

The cyclomatic number has several related interpretations: it counts the number of 
independent cycles in the graph (hence the name) and therefore it is the first Betti 
number of T (hence the notation j3). It also counts the minimal number of edges that 
need to be removed from T to turn it into a tree. Correspondingly, /5p = if and only 
if T is a tree. 
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There is another simple but useful observation relating the cycles on the graph and the 
number of zeros: if the eigenfunction of the graph does not vanish on the vertices of the 
graph, the number of zeros on any cycle of the graph is even. Indeed, an eigenfunction 
of a second order operator can only have simple zeros, thus at every zero / changes 
sign. On a cycle there must be an even number of sign changes. 

As mentioned earlier, we will be interested in the number of zeros and nodal domains 
of the eigenfunctions of operators (11. ip and (jl.2p on graphs. According to the well 
known ODE theorem by Sturm [32| |33| [31], the zeros of the n-th eigenfunction of the 
operator of type ( )1.2p on an interval divide the interval into n nodal domains. By 
contrast, in the corresponding question in IR 2 only an upper bound is possible, given by 
the Courant's nodal line theorem [5], v n < n. In a series of papers |14| [T5| [T%j IT?] 20J, 
it was established that a generic eigenfunction of the quantum graph satisfies both an 
upper and a lower bound. Namely, let A n be a simple eigenvalue of — -pr + V(x) on a 
graph T and its eigenfunction /„ be non-zero at all vertices of T. Then the number of 
the nodal domains of f n satisfies 

(1.9) n - p T < u n < n. 
Similarly, for the number of zeros we have 

(1.10) n- 1 < n n <n- l + /3 r . 

Note that the upper bound in (11.101) follows from the upper bound in (jl.9p and inequality 
(II. 7p . The lower bound in (ll.lOp requires an independent proof which is given in |35| . 
An interesting feature of quantum graphs is that, unlike the M. d case, the upper bound 
v n < n is in general not valid for degenerate eigenvalues. 

In the present paper we combine these a priori bounds with scattering properties of 
a certain family of graphs to derive formulas for the nodal counts fi n and u n . 

1.3. Quantum evolution map. When the potential V(x) is equal to zero, the eigen- 
value equation 

dx 1 

has, on each edge, a solution that is a linear combination of the two exponents e ±lfcr if 
k ^ 0. We will write it in the form 

(1.12) fe{x e ,v) = a^v exp(-ikx e>v ) + a™' exp(ikx e , v ), 

where the variable x e>v measures the distance from the vertex v of the edge e. The 
coefficient a e l ™ is the incoming amplitude on the edge e (with respect to the vertex v ) 
and a e °"* is correspondingly the outgoing amplitude. However the same function can be 
expressed using the coordinate 2 e u as 

(1.13) f e {%e,u) = a e m u exp(-i/cx eiU ) + ag£ exp(ikx e>u ). 

Since these two expressions should define the same function and since the two coordi- 
nates are connected, through the identity x e>v + x e , u = L e , we arrive to the following 
relations 



(1.14) a% = e lkL ^ a:% = 

Fixing a vertex v of degree d v and using f l 1 . 1 2 j) to describe the solution on the edges 
£ v adjacent to v we obtain from ( 11. 3p and (JL4]) d v equations on the 2d v variables a™ f 
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and a*™. These equations can be rearranged as 

(1.15) a° ut = ot"\k)a™, 

where a° ut and a* n are the vectors of the corresponding coefficients and o~( v \k) is a 
d v x d v unitary matrix. The matrix a^°'(k) is called the vertex-scattering matrix, it 
depends on k for values of a v other than or oo and its entries have been calculated 
in [36] . 

Collect all coefficients a]" into a vector a of size 2\E\ and define the matrix J acting 
on a by requiring that it swaps around a*™ and a*" for all e = (u, v). Then, collecting 
equations fll . 15j) into one system and using connection (11.141) and the matrix J to rewrite 
everything in terms of a we have 

Je~ [kL a = E(fc)a. 

Here all matrices have the dimension equal to double the number of edges, 2\S\. The 
matrix L is the diagonal matrix of edge lengths, each length appearing twice and £ is 
the block-diagonalizable matrix with individual trW as blocks, namely 

Zj (ei,?;i),(e2,t>2) ~~ °v 1 ,V2 <J ei ,ei " 

Noting that J -1 = J, the condition on a can be rewritten as 

(1.16) a = e ikL JZ(k)a, 

The unitary matrix U(k) := e lfcL JS(fc) is variously called the bond scattering matrix 
[36J or the quantum evolution map [2J. The matrix JS(fc) describes the scattering of 
the waves on the vertices of the graph and e lhL gives the phase shift acquired by the 
waves while traveling along the edges. The quantum evolution map can be used to 
compute the non-zero eigenvalues of the graph through the equation 

(1.17) det (I — U(k)) = 0. 

We stress that U (k) is not a scattering matrix in the conventional sense, since the graph 
is not open. Turning graph into a scattering system is the subject of the next section. 



2. Attaching infinite leads to the graph 

A quantum graph may be turned into a scattering system by attaching any number of 
infinite leads to some or all of its vertices. This idea was already discussed in [36 | [24^ , 137]. 
We repeat it here and further investigate the analytic and spectral properties of the 
graph's scattering matrix, that would enable the connection to the nodal count. 

Let T = (V, £) be a quantum graph. We choose some M < |V| out of its vertices 
and attach to each of them an infinite lead. We call these M vertices, the marked 
vertices, and supply them with the same vertex conditions as they had in T. Namely, 
each marked vertex v retains its 5-type condition with the same parameter a v (recall 
(II. 4p ). We denote the extended graph that contains the leads by T and investigate its 
generalized eigenfunctions. 

The solution of the eigenvalue equation, (II. lip , on a lead / which is attached to the 
vertex v, can be written in the form 

(2.1) fi(xi,v) = ci n v exp(-ikx i:V ) + cffexpfax^). 

The variable xi v G [0, oo) measures the distance from the vertex v along the lead I 
and the coefficients c/", c™ ( are the incoming and outgoing amplitudes on the lead I 



NODAL COUNT ON A FAMILY OF QUANTUM GRAPHS 



7 



(compare with (I1.12p ). We use the notation c out , c in for the vectors of the corresponding 
coefficients and follow the derivation that led to (I1.16P in order to obtain 



(TjOUt \ 
° E )=e ikL °J o m 

All the matrices above are square matrices of dimension 2 |£|+M. There are two differ- 
ences from equation ( I1.16p . First, in the matrix L a each lead is represented by a single 
zero on the diagonal, in contrast to the positive lengths of the graph edges, appearing 
twice each. The matrix J a swaps around the coefficients a corresponding to opposite 
directions on internal edges, but acts as an identity on the leads. These differences 
arise because for an infinite lead we do not have two representations fll . 12j) and f ll . 13[) 
and therefore no connection formulas fll,14j) allowing to eliminate outgoing coefficients. 
Writing the matrix e lkL JE(fc) in blocks corresponding to the edge coefficients and lead 
coefficients results in 



(2.3) 



c out \ _ ( R{k) T (k) 
T,{k) U{k) 



where the dimensions of the matrices R, T Q , Tj and U are M x M, M x 2\S\, 2\S | x M 
and 2\S\ x 2\£\ correspondingly. We stress that the matrix U (k) describes the evolution 
of the waves inside the compact graph and has eigenvalues that can now lie inside the 
unit circle due to the "leaking" of the waves into the leads. 

Equation (12. 3p can be used to define a unitary scattering matrix S such that c out = 
S c in , as described in the following theorem. 

Theorem 2.1. Let 

(2.4) ( C 7 )=«(£). ^ere Q = ( £ £ 

is a unitary matrix with the blocks R, T , Ti and U of sizes M x M, M x 2\S\, 2\S\ x M 
and 2\S\ x 2\S\ correspondingly. For every choice of c m G C M , consider relation ( 12.4ft 
as a set of linear equations in the variables c out G C M and a G C 2 ' £ '. Then 

(1) There exists at least one 2\S\ x M matrix C such that 

ii 



(2.5) (I - V) C = Ti 

(2) Let 

(2.6) S = R + T C. 

Then S is a unitary matrix independent of the particular choice of C in equa- 
tion (123]) . 

(3) The solutions of (12. 3p are given by 

(2.7) c out = Sc 111 , 

(2.8) aG {Cc in + Ker (I- U)}. 

In particular, c out is defined uniquely by c in . 

The proof of the theorem distinguishes between the case of a trivial Ker(I — U) and 
the case of singular I — U. The following lemma makes the treatment of the latter case 
easier. 
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Lemma 2.2. Let Q be as in theorem \2.1[ Then the following hold: 



(2.9) 
(2.10) 



Range Tj C Range (I — U), 
Ker(I- U) C KerT 



Proof. Since in a finite-dimensional space Range A = (Ker A*) 
alent to 



equation (12. 9p is equiv- 



(KerT*)^ Q (Ker (I - [/*) 
which is in turn equivalent to 

Ker (I -IT) C Ker T*. 
Let v G Ker (I — U*). Using the unitarity of Q* we get 



(2.11) 



R* T* 
T* U* 



T*v 
U*v 



v 



Equating the left-hand side to the right-hand side of the equation above we get T*v = 0. 
Equation (I2.10p is proved in a similar manner by replacing Q* with Q in the above. □ 

Proof of theorem\2Ji Case 1: det(I -U)j= 0. 

To show part ([T|) we simply set C = (I — U)~ x T{. Furthermore, equation (12.31) has a 
unique solution, given by 

-l 



(2.12) 
(2.13) 



I — U) 



tout 



(R + T (l-U) l T^j 



which proves part ([3]). 

The unique definition of C guarantees the uniqueness of S = R + T C. To finish the 
proof of part fl2J) we use the unitarity of Q, which provides the identities 

R*R + T*Ti = T*T Q + U*U = I 

(2.14) R*T + T*U = T*R + U*Ti = 0. 

From here we get 

S*S = R*R + R*T Q C + C*T*R + C*T*T C 

= I - T*Ti - T*UC - C*U*Ti + C*(l- U*u) C. 



Expanding, factorizing and using the definition of C in the form Tj + UC = C, we 
arrive to 

S*S = I + C*C - {T* + C*U*)(Ti + UC) = I 
Case 2: det(I - U) = 0. 

Existence of a solution C to the equation (1—U)C = T\ is guaranteed by equation (12. 9p 
of Lemma 12.21 

The columns of C are defined up to addition of arbitrary vectors from Ker (I — U). 
However, Lemma 12. 2\ equation ( I2.10p implies that these vectors are in the null-space 
of T , therefore the product T Q C has unique value independent of the particular choice 
of the solution C. The proof of the unitarity of S has already been given in case 1 and 
did not rely on the invertibility of I — U. This proves part (J2j). 
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The last 2\E\ equations of ( 12. 3ft are (I — U)a = TiC %n . From (12. 5ft . all solutions a of 
this equation are given by Cc in + Ker(I — U). On the other hand, the first M equations 
of ( 12. 3p are c out = Rc m + T a and substituting the already obtained expression for a 
and using ( 12 . 1 0[) we finally arrive to c out = (R + T C)c m . This finished the proof of the 
theorem. □ 

We would like to study the unitary scattering matrices, S(k) as a one-parameter 
family in k G R. The matrix Q(k) is a meromorphic function of k in the entire complex 
plane (29]. For all k values which satisfy det(I — U) 0, S(k) is given explicitly by 

~ \ -i 



(2.15) S(k) = R + T (l - U 



T, 



i ■ 



and S(k) is therefore also a meromorphic function in k at these values. The significance 
of the values of k for which det(I — U) — is explained in the following lemma. 

Lemma 2.3. Let T* be the quantum graph obtained from the original compact quantum 
graph, T, by imposing the condition f(v) = at all of its M marked vertices, in addition 
to the conditions already imposed there. Then the spectrum a (T*) coincides with the 
set 

(2.16) A={*; 2 gM det(l-C/)=0 

Proof. We mention that imposition of the additional vertex conditions makes the prob- 
lem overdetermined. In most circumstances the set A will be empty. The operator H 
is still symmetric but no longer self-adjoint, because its domain is too narrow. 

Denote by V the graph with the leads attached. Let k 2 G a (r*) and let / be the 
corresponding eigenfunction on T*. Then / can be extended to the leads by zero. It will 
still satisfy the vertex conditions of T and will therefore satisfy (12. 3p with c m = c out = 
and a ^ 0. The last 2\£\ equations of ( 12. 3 j) imply det(I — U) = 0. 

In the other direction, let k 2 G A. Choose a G Ker(I — U). We see that equation 
(12. 3p is satisfied with the chosen a and with c m = c out = 0. These coefficients describe 
a function on T which vanishes completely on the leads and therefore its restriction to 
T satisfies Dirichlet boundary conditions on the marked vertices by continuity. This 
implies k 2 G a (T*). □ 

Corollary 2.4. The set A is discrete. 

Proof. The corollary is immediate since A C o~(T), which is discrete. □ 
Lemma 2.5. S (k) is a meromorphic function which is analytic on the real line. 



Proof. The blocks of the matrix in equation (12.2ft are meromorphic (see [29], 

Theorem 2.1 and the discussion following it), therefore all the blocks of the matrix Q 
are meromorphic on the entire complex plane. Since the set A on which the matrix I—U 
is singular is discrete, equation (I2.15P defines a meromorphic function. To show that 
S(k) in fact does not have singularities on the real line, we observe, that, for k G R\ A 
we have shown that S(k) defined by (I2.15P is unitary. Therefore S(k) remains bounded 
as we approach the "bad" set A and the singularities are removable. Theorem 12. 1 1 gives 
a prescription for computing the correct value of S(k) for k G A. □ 

We now examine the ^-dependence of the eigenvalues of S(k). To avoid technical 
difficulties we restrict our attention to the case when only a = (Neumann) or a = oo 
(Dirichlet) are allowed as coefficients of the <5-type vertex conditions, equation (jl.4p . In 
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this case the matrix described in section [L~3l is independent of k making calculations 
easier. The general case can be treated using methods of [38j, however we will not need 
it for applications. 

Lemma 2.6. Let every vertex of the graph V have either Neumann or Dirichlet con- 
dition imposed on it. Then the eigenvalues of S(k) move counterclockwise on the unit 
circle, as k increases. 

Proof. Let e lLp be an eigenvalue of S with the normalized eigenvector v. Differentiating 
the normalization condition v*v = 1 with respect to k we get 

(2.17) v*v = 0. 

Now we take the derivative of Sv = e tip v with respect to k to get 

d 



S)v + Sv = iCpe^v + e^v. 
dk / 

We multiply the above equation on the left with v*S* = e~ ltp v* and use v*v = 1 and 
equation fl2. 1T[) to obtain: 

v b —b ) v = up. 
\dk J 

Thus we need to show that \S*4rS is positive definite. Comparing equations (I2.2p and 
(I2.3p and using that X is /^-independent, we obtain that 

R{k) = R(0), T (k) = T o (0), Ti(k) = e ifcL T,(0), U(k) = e ikL U(0). 
Differentiating the latter two matrices with respect to k produces 

U = iLU, fi = iLTi. 
We can now differentiate equation (I2.5P to obtain 

(2.18) (1 - U)C = fi + UC = iL{Ti + UC) = iLC, 

where we used (12. 5p again in the final step. 
For the matrix in question we now obtain 

b*4rS = (R* + C*T*) T C = -T*UC + C*(l- U*u) C, 
dk \ / 

where equations (I2.14p have been used in the second step. Using —T* = C*(U* — I) 
which is a conjugate of (|2.5|) . we obtain 

S*-^b = C* (p*U — U + I— U*U^j C. 
Using (I2.18P this simplifies to 

S* —b = iC*LC. 
dk 

Since L is diagonal with positive entries we conclude that jS^^S* is positive definite. □ 

We end this section by stating a result known as the inside-outside duality, which 
relates the spectrum of the compact graph, T, to the eigenvalues of its scattering matrix, 
b(k). This is a well known result, mentioned already in We bring it here with a 
small modification, related to the already mentioned set, A. 

Proposition 2.7. The spectrum ofT is A U {k \ det (I — S) = 0}. 
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Proof. We remind the reader that when a lead is attached to a (marked) vertex, the 
new vertex conditions are also of 5- type with the same value of the constant a v . The 
conditions at the vertices that are not marked remain unchanged. 

Let k be such that det (I — S) = 0. Let c be the corresponding eigenvector, Sc = c. 
Letting c m = ewe find c out = c and a according to theorem 12.11 The corresponding 
generalized eigenfunction / satisfies correct vertex conditions at all the non-marked 
vertices. It is also continuous at the marked vertices and satisfies 

(2.19) E £(„ )+E ^ W=aJW , 

where S v is the set of the finite edges incident to v and is the set of the infinite leads 
attached to it. Referring to (12. ip we notice that c*™ = c™* implies that the derivative of 
/ on the lead is zero. Therefore equation ( I2.19P reduces to the corresponding equation 
on the compact graph. Thus the restriction of / to the compact graph satisfies vertex 
conditions at all vertices and k 2 is an eigenvalue of T. Inclusion A C a(T) has already 
been shown in lemma |2~3"1 

Conversely, let k 2 be an eigenvalue of the compact graph T and let / be the corre- 
sponding eigenfunction. Then / can be continued onto the leads as f(v) cos(kx), where 
f(v) is the value of / at the vertex v where the lead is attached to the graph. Compar- 
ing to (12. ip we see that c\ n v = c° u J = f(v)/2. Therefore the resulting extended function 
is characterized by vectors a, c m and c out such that c m = c out = Sc m . If the function 
/ was non-zero on at least one of the marked vertices, c m ^ is a valid eigenvector of 
S(k) with eigenvalue 1. If / is zero on all marked vertices, k e A by lemma [2731 □ 



3. Applications to the nodal domains count 

3.1. Application for a single lead case. We wish to study the nodal count sequence 
of a certain graph V by attaching a single lead to one of its vertices. Let S(k) = e ilf ^ k ' 
be the corresponding one dimensional scattering matrix. For each real k > there 
exists a generalized eigenfunction, f(k; x), of the Laplacian with eigenvalue k 2 on the 
extended graph, T, as proved in theorem 12. II In addition, up to a multiplicative factor, 
this function is uniquely determined on the lead, where it equals 

fiead(k; x) = c m exp(-ikx) + c out exp(ikx) = c m [exp(— ikx) + exp (i (<p(k) + kx))] 



c exp i 



m Q0S ( kx+ m 



2 J \ 2 J 



The positions of the nodal points of this function on the lead are therefore uniquely 
defined for every real k > and given by 



(3.1) D lead (k) = \x>0 



09 (k) 71 7T „ 

We exploit this by treating A; as a continuous parameter and inspecting the change in 
the positions of the nodal points as k increases. Let x = — + ^ + ^n x be the 
position of a certain nodal point on the lead at some value k, i.e., x G D\ ead (k). The 
direction of movement of this nodal point is given by 
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d / tp (k) 71 71 \ klfi'(k) — (p(k) + 7T + 27TTl x 

dk V 2k~ + 2k + k Ux ) = 2k 2 

where for the last inequality we need to assume that all the vertex conditions of V are 
either of Dirichlet or of Neumann type in order to use the conclusion of lemma 12.61 , 
if' (k) > 0. From (13. 2p we learn that all nodal points on the lead move towards the 
graph, as k increases. 

The event of a nodal point arriving to the graph from the lead occurs at values k 
for which G Di ea d (k) and will be called an entrance event. We may use (13. ip to 
characterize these events in terms of the scattering matrix: 

(3.3) G D lead (k) mod 27r {f{k)) = tt. 

After such an event occurs, the nodal point from the lead enters the regime of the graph 
T and may change the total number of the nodal points of f(k; x) within V. 
Another significant type of events is described by 



(3.4) mod 27r {(f{k)) = 0. 

Proposition 12.71 shows that such an event happens at a spectral point of T and at this 
event, the restriction of f(k; x) on V equals the corresponding eigenfunction of V . These 
events form the whole spectrum of V if and only if A = 0. This is indeed the case if we 
choose to attach the lead to a position where none of the graph's eigenfunctions vanish 
(lemma [273]) . In addition, we will assume in the following discussion that V has a simple 
spectrum. This is needed for the unique definition of the nodal count sequences and is 
shown in [3H] to be the generic case for quantum graphs. 

The two types of events described by (13.31) and (13. 4p interlace, as we know from 
lemma 12.61 (compare also with theorem lA.ip . We may investigate the nodal count 
of tree graphs by merely considering these two types of events and their interlacing 
property. We count the number of nodal points within V only at the spectral points 
to obtain the sequence {/i n }- Between each two spectral points we have an entrance 
event, during which the number of nodal points within V increases by one, as a single 
nodal point has entered from the lead into the regime of V. This interlacing between 
the increments of the number of nodal points and its sampling gives /i n = n — 1 and 
v n = + 1 = n. 

The above conclusion is indeed true for tree graphs under certain assumptions (see [HI 
[T51[T7]). However, when graphs with cycles are considered, there are other interesting 
phenomena to take into account: 

(1) In the paragraph above it was taken for granted that at an entrance event the 
nodal points count increases by one. This is indeed so if the nodal point which 
arrives from the lead enters exactly one of the edges of T without interacting 
with other nodal points which already exist on the graph. However, when the 
lead is attached to a cycle of the graph, the generic behavior is either a split or a 
merge. Assume for simplicity that the attachment vertex has degree 3, counting 
the lead. A split event happens when a nodal point from the lead splits into two 
nodal points that proceed the two internal edges. This will increase the number 
of nodal points on T correspondingly. In a merge event the entering nodal point 
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merges with another nodal point coming along one of the internal edges. The 
resulting nodal point proceeds along the other internal edge. The number of 
nodal points on T will not change during such an event. If a lead is attached to 
a vertex of higher degree the variety of scenarios can be larger. 
(2) Another type of events that were not considered are ones in which a nodal point 
travels on the graph and reaches a vertex which is not connected to the lead. 
When this vertex belongs to a cycle, the generic behavior would be a split or 
a merge event and would correspondingly increase or decrease the number of 
nodal points present inside the graph. 

These complications are dealt with in section H~Tj where we use the single lead approach 
to derive a nodal count formula for graphs which contain a single cycle. In the following 
section, 13.2} we consider a modification of this method — we attach two leads to a graph 
and use the corresponding scattering matrix to express nodal count related quantities 
of the graph. Later, in section 13. 'S\ we show how the two leads approach yields an exact 
nodal count formula for a specific graph. 

3.2. Sign-weighted counting function. The number of nodal points on a certain 
edge e = (u,v) at an eigenvalue k n is given by 



where [^J stands for the largest integer which is smaller than x, and f n is the cor- 
responding eigenfunction [18J. We infer that the relative sign of the eigenfunction at 
two chosen points is of particular interest when counting nodal domains. While the 
most natural candidates for the two points are end-points of an edge, the results of this 
section apply to any two points on a graph. Denoting these points x\ and X2, we are 
interested in the sign of the product f n {xi)fn{x2), where /„ is the n-th eigenfunction 
of the graph. We define the sign-weighted counting function N XlyX2 (k) as 



(3.6) N XuX2 (k) = #{k n < k : f n (x 1 )f n (x 2 ) > 0} - #{k n < k : f n {x{)f n {x 2 ) < 0}. 



Using the scattering matrix formalism allows us to obtain the following elegant formula 

Theorem 3.1. Let V be a graph with Neumann or Dirichlet vertex conditions and x\ 
and x 2 be points on the graph such that no eigenfunction turns to zero at X\ or x 2 - 
Denote by S(k) the 2x2 scattering matrix obtained by attaching leads to the points x\ 
and X2 ■ Let I e be the matrix 



where e 2l ^ k ' = det(S(k)) and a suitable continuous branch of the argument is chosen. 
The convergence is pointwise everywhere except atkE cr(Y) (where N Xl , X2 is discontin- 
uous). 

Proof. First of all, we observe that the scattering matrix of a graph with only Neumann 
or Dirichlet conditions is complex symmetric: Sj^ = Skj- This can be verified explicitly 
by using representation (I2.5|) - (l2.6p . together with (I2.2D and the fact that the matrix £(&) 
is real and symmetric under the specified conditions. 



(3.5) 




(3.7) 





lim arg det(I e - S(k)) - (p(k) , 
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For a moment, consider that only one lead is connected, to the point x\. Then the 
events that change the sign of f{x\) and f(x 2 ), where / is the one- lead scattering 
solution, are the values of k such that (A) a zero comes into the vertex x\ from the 
lead ("Dirichlet events") or (B) a zero crosses the point x 2 . The former are easy to 
characterize: they interlace with the events when a "Neumann point" comes into the 
vertex X\, which happen precisely at cr(T), as discussed in section |3~T1 

Denote by k the value of k when a zero crosses the point x 2 where the lead is not 
attached. Now consider the scattering system when both leads are attached, at points 
Xi and x 2 . At k = k the one- lead scattering solution f(x) can be continued to the 
second lead by setting it to vanish on the entire lead. This would create a valid two- 
lead solution with c 2 n = c ^ 1 = 0. By inspecting (12. 7p we conclude that the vector 
(1, 0) T is therefore an eigenvector of the two-lead scattering matrix S(k). This happens 
whenever the matrix S(k) is diagonal. We conclude that the events of type (B) happen 
in a one-lead scattering scenario precisely when the two-lead scattering matrix satisfies 
S(k\ 2 = S{k) 2>1 = 0. 

Introducing the notation, 

(3.9) ((k) = det(J - S(k)) 

(3.10) r{k) = S(k) li2 = S(k) 2>1 , 

we can summarize the earlier discussion as follows. The eigenvalues of the graph are 
given by the zeros of ((k) (see proposition 12. 7p . and the relative sign of the n-th eigen- 
function, sign[/ n (a;i) f n {x 2 )}, is equal to the parity of the total number of zeros of £ and 
r that are strictly less than k n . Note that the condition that no eigenfunction is zero 
on Xi or x 2 implies that the set A in proposition 12.71 is empty and that the zeros of the 
functions ( and r are distinct. 

Applying complex conjugation to ((k) we obtain 

(3.11) C(kY = det(S*)((k) = C(fc)/det(S). 

Similarly, using the explicit formula for the inverse of a 2 x 2 matrix together with the 
unitarity of 5*, we obtain for r(k) 

(3.12) r(k)* = - det(S*)r(k) = -r(Jfe)/ det(S). 

These relations allow us to represent ((k) = r{k)e l( ^ <yk ^ and r(k) = -t^e 1 ^ 1 *' , when 
recalling that e 2im = det(S(k)). 
We now evaluate 

det(I £ - S) = det(I - S) - e{S x , 2 + S 2A ) - e 2 

= ((k) - er{k) -t 2 = (r(k) - iet(k))e im - e 2 , 

and, therefore, 

(3.13) det(I e - S)e~ im = {r{k) - iet(k)) + o(e). 

It is now clear that, when r(k) ^ (i.e. when k (jL o~(T)), the e — > limit of the above 
ratio is a non-zero real number and therefore its argument is an integer multiple of n. 

To evaluate this integer we focus on the values of k when z(k) = r(k) —iet(k) crosses 
the line Re(z) = 0. When the crossing is in the counter-clockwise direction, the integer 
above increases, and otherwise it decreases. The counter-clockwise versus clockwise 
direction of the crossing is decided exclusively by the sign of the ratio r{k n — 0)/t(k n — 0), 
which coincides with the parity of the total number of zeros of the two functions. This, 
in turn, has been shown to coincide with the relative sign of the eigenfunction. □ 
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Remark 3.2. It is interesting to compare the above formula for the sign-weighted count- 
ing function with the corresponding formula for the more commonly used spectral 
counting function, 

(3.14) N(k) = #{k n < k}. 

Under conditions of theorem 13. II the counting function N(k) can be represented as 

(3.15) N(k) = - flimargdet(I - el - S) - <f>(k) 

Combining the two we can obtain a counting function that counts only the eigenvalues 
whose eigenfunctions have differing signs at X\ and x 2 , 

(3-16) N~ uX2 (k) := #{k n < k : f n (x 1 )f n (x 2 ) < 0} = ±- limarg d< * (1 " cI " S) 

2n e->o det(I — eo\ — b) 

3.3. Using two leads to derive an exact nodal count formula. In the current 
section we derive a nodal count formula for the graph, T, given in figure [37T1 This graph 

2b 

vN 

n r > d 




2c 

FIGURE 3.1. The quantum graph whose nodal count we compute. The 
lengths of the edges and the vertex conditions are indicated. D stands 
for Dirichlet vertex conditions and tV for Neumann ones. 



is a member of an isospectral pair, as described in [40J. The isospectral twin of this 
graph is the graph shown in figure | 



2a 





FIGURE 3.2. The tree graph which is isospectral to the graph in figure 
13.11 The lengths of the edges and the vertex conditions are indicated. 

Examining the topology of each of the graphs according to (jl.9p tells us that the tree 
graph has the nodal count v n = n and the nodal count of the graph T, which has a 
cycle, obey the bounds n — 1 < v n < n. It was claimed in |41| that the nodal count of 
T is 

(3.17) v n = n- X - L ^" J . 

This formula was not proved there, but rather a numerical justification was given. We 
present here a proof for the following theorem 

Theorem 3.3. Let a, b, c be positive real numbers such that | ^ Q and ^ ^ Q. Let T 
be the graph described in fiaure \3.1[ Then the nodal points count sequence ofT is 



(3.18) fi n = n + mod2 



b + c 

; n 

a + b + c 
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and the nodal domains count sequence of T is: 



(3.19) 




n 



n — 1 + mod 2 ( |_ 



a+6+c 



b+c 



n 



» < L^J + 1 

J) n >lWc\+ 1 



Remark 3.4. Note that for eigenvalues high enough in the spectrum, the exact nodal 
count (I3.19P coincides with the previously known numeric result, (I3.17p . 

The method of proof of the formulas (I3.18p . (13 . 1 9[) involves attaching leads to the 
graph and keeping track of the nodal points dynamics with respect to an increment of 
the spectral parameter. This specific example both presents the ability to derive an 
exact formula and also demonstrates the technical complications that may arise while 
using this method. 

3.3.1. A brief outline of the proof. Define a graph with two vertices connected by two 
edges of lengths 2b and 2c. This graph is actually a cycle of length 2b + 2c. Connect 
two leads to the vertices of this graph to obtain the graph in figure 13.31 Denote this 
graph by T and notice that T as a metric graph is a subgraph of T. The graph T has 
a symmetry of reflection along an axis which passes through the middle of the graph. 
We will exploit this symmetry in the next section. 



Figure 3.3. The graph T and its symmetry axis (dotted). The values 
of its generalized eigenfunction are specified on the leads. 

The Laplacian on V possesses a continuous spectrum and each generalized eigenvalue, 
k 2 , has a two-dimensional generalized eigenspace, characterized by c m G C 2 (theorem 
12. ip . We will describe a one parameter {k G R) family of generalized eigenf unctions on 
T, f(k; x). We thus consider f(k; x) as a function on V which changes continuously 
with k - this will be emphasized by the notation f(k; •). This ^-dependent function 
would be chosen such that its restriction on the subgraph T at k 2 G a (V) equals the 
corresponding eigenfunction of T. The strategy of the proof is to keep track of the 
number of nodal points of f(k; ■) as it changes with k and to sample this number 
at k 2 G cr(r). We will notice that the nodal points travel continuously from infinity 
towards the cycle and we will characterize the dynamics of the nodal points which 
enter the cycle. This will allow us to find the change in the number of nodal points 
during such entrance events. We will then calculate the number of eigenvalues which 
occur between two consequent entrance events and will combine all those observations 
to deduce the nodal count formulas (I3.18P and (I3.19p . 

3.3.2. Towards a proof of theorem \3.3[ Let F be the graph that is described in section 
13. 3. H and appears in figure I3~3"l A generalized eigenfunction of T with eigenvalue k 2 , on 
the ] th lead, is given by 

fj(k; Xj) = cj 1 exp(-ikxj) + c° ut exp(ikxj), 



cf* exp (— ikxi) + c° ut exp (ikxi) 



CD 



eg 1 exp (— ikx2) + c% ut exp (ifca^) 
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where j = 1,2 and the coefficients 



V c ? J v c °2 ut 

are related by 

(3.20) c out = Sik)^. 

The graph T obeys a symmetry of reflection along a vertical axis which passes through 
the center of the graph. This reflection symmetry exchanges the two leads of T and it 
implies that its scattering matrix, S(k), commutes with the matrix 

1 

1 



.out 



a 



This, together with the unitarity of S(k) (theorem 12. ip allows us to write it in the form 

cos(7(A;))e^ fc ) sin (7 (k)) e*M k )+ K / a ) 
sin (7 (k)) e^( fc ) + "/ 2 ) cos (7 (Jfe)) e 1 ^ 



(3.21) S(k) 



The exact form of S(k) (expressed in terms of the edge lengths parameters a, b, c) can 
be calculated using fl22J , (J2~3D and (j2~gj) . 

Following the approach described in section I3~T1 we treat A; as a continuous parameter 
and choose c in to vary continuously with k. Namely, we choose a certain continuous 
vector function c m : (0, 00) — > C 2 . Relation ( 13 .20 j) yields the continuous function 
c out : (0, 00) — > C 2 and both c m (k) and c out (k) determine f(k; •), a function on T that 
changes continuously with k. We next describe a specific choice of c m (k) that yields a 
function f(k; x) with the following properties which are convenient for our proof. 

Property 3.5. The values of the function on the leads are real, i.e., 



f j (k;x j ) = f j (k;x j ) for j = 1,2. 
Property 3.6. Denote the zeros of the function and of its derivative on the leads by 

Dj(k) := {,->() ./•(/,•:,;,) } 

f d 

Nj(k) := \ Xj >0 



fj(k; xj) = 



They obey 

D 1 (k) = N 2 (k) and D 2 (k) = N t (k) . 
The usefulness of these properties is made transparent in the following proposition. 

Proposition 3.7. Let a,b,c be positive real numbers, such that - ^ Q. Let T and T 

the graphs defined above (with the edge lengths parameters a,b,c). 

(1) For each fceR properties \ 3.5\ and \3.6i define a function f(k; •) on T which is 
unique up to a multiplication by a scalar and a reflection along a vertical axis 
which passes through the middle ofY. 

(2) The above function, f(k; ■), can be chosen to be continuous in k. 

(3) If k 2 G o~(T), the restriction of the function f(k; •) to the graph on T coincides 
with the eigenfunction of T up to reflection. 

The following lemma will aid us in proving the uniqueness of f(k; •). 
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Lemma 3.8. Let b, c be positive real numbers, such that - ^ Q. Then the set 



A = Ue 



det I - U = 



that was defined in (12. 1 6[) . lemma WT^. is an empty set. 

Proof. Lemma [2.31 tells us that A = o (T*), where T* is a cycle graph with additional 
Dirichlet conditions imposed on its two vertices (figure [33]). 

2b 




2c 

FIGURE 3.4. The graph P", which obeys A = cr(T*). 

Assume that k is in the spectrum of T*. The corresponding eigenfunction should 
then be of the form sin (kx) on each of the edges (up to a multiplicative factor). The 
Dirichlet boundary conditions imply that sin (2bk) = sin (2ck) = and therefore b and 
c both belong to the set ^N. This means that ' 6 Q and contradicts the assumption. 

□ 

Proof, of proposition 13.71 

Let k G ffi. Let f(k; •) be a generalized eigenfunction of T which obeys the properties 
13.51 and 13.61 From property 13.51 we conclude that 



c f = df for j = 1,2. 



Thus, for a suitable Cj and <x,-, 



in / • * 

Cj = Cj exp yiaj i 



(3.22) c° ut = Cj exp (-i aj ) . 

We plug this in the expression for the values of / on the leads 

fj(k; Xj) = cj 1 exp(— ikxj) + c° u * exp(ikxj) = 2cj cos(<x,- — kxj), 

and obtain 

Dj (k) = {Xj > \fj(k; Xj) = 0} 



Xj > 



Property 13.61 now translates to 
(3.23) 



Xj > 



Xj > 



_d_ 

dxj 



fj(k, Xj 



«2 



7T 



7T 



z 



We use f l3.2ip .f l3T22l) and (I3.23f) and plug them in (13.201) to get equations on Cj,aj. 
There are two possible solutions, which describe two functions that are the same up 
to a reflection along a vertical axis which passes through the middle of T. One of the 
solutions reads 
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(3.24) c\ n (k) = ^exp(-^ (fc) 



2 1 \ 2 

(3.25) c?(k) = 



cos(7(fc))-l ovT ^ ( _■ y(fc)-?r 
2sin(7(fc)) 



exp f-i^hA 7 (k) i 2ttZ 
7 (Jfc) £ 2vrZ 

and the corresponding function is given on the leads by 



xi) = cos ( — V kxi 



,£ Sjlr-s(^ + fcx 2 ) 7 (*)*2ttZ 
7(A;)g2ttZ 



Note that f2{k; ■) is continuous in k. In addition, ■) and f 2 {k] •) that are given 

above can be multiplied by any ^-continuous scalar function to yield an appropriate 
solution which is also continuous in k. This proves that f(k; •) is uniquely defined on 
the leads and also fc-continuous there. It is left to show the same for the values of 
f(k; •) on the cycle. Theorem 12.11 implies that f(k; •) may have multiple values on the 
cycle only for k 2 £ A. However, since A = (lemma l3.8p . this cannot happen and 
f(k; ■) is uniquely defined on the cycle. In addition, the values of f(k; •) on the cycle 
are determined by equation ( 12.121) . which shows that these values are continuous in k, 
due to the reversibility of I — U and the ^-continuity of c in . 

We start proving part [3] of the proposition by assuming that k 2 £ cr(T). We have 
that there exists a real eigenfunction with eigenvalue k 2 on T. We fix a function g(k; •) 
on T to equal this eigenfunction when restricted on F. Then the values of this function, 
g(k; •), can be uniquely continued so that it is defined on the whole of T. It is easy to 
verify that the obtained function obeys properties 13.51 and 13.61 and we conclude from 
the proof of part [TJof the proposition that it is equal to f(k; •) up to a multiplication 
by a scalar or a reflection. 

□ 

Proposition 13 . 71 shows that there are only two /^-continuous functions, f(k; x), which 
obey the properties 13.51 and 13.61 We call such a function a real contra-phasal solution, 
due to the properties that it has. These functions will be used to prove theorem 13.31 
We carry on by stating a few lemmas which describe the dynamical properties of the 
nodal points of such a real contra-phasal solution. 

Lemma 3.9. The nodal points of a real contra-phasal solution move on the leads to- 
wards the cycle as k increases. 

Proof. While proving proposition 13.71 we have showed that one of the real contra-phasal 
solutions has the following values on the leads 

fi(k; xi) = cos ( — ^- + kx\ 



(3.26) f 2 {k- x 2 ) 



' £ S^rCos(^ + fcx 2 ) 7 (*)£2*Z 
) 7(A;)£2ttZ 
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The positions of its nodal points on the leads are therefore given by 
Di(Jfe) = |xi>0 



(3.27) D 2 (k) = <x 2 >0 



09 (k) 71 71 

Let x (k) = — + ^ + \n x be the position of a certain nodal point on the first 
lead at the value k, i.e., x (k) G D\ (k). The direction in which this nodal point travel 
on the first lead is given by 

, ... d / cp (k) 7i 7i \ kip'(k) — <p(k) + tx + 27™ 

x' (k) = — -^-^r 1 + — + T n 1 ' 



(3.28) = -U^+x 



dk V 2k ' 2k ' k'" x J 2k 2 

1 
k 

A simple calculation based on (13.21)) gives 

det S (k) = exp (i2ip (k)) . 

Denoting the eigenvalues of S (k) by exp (iipi (k)) , exp [i , p 2 {k)), we have that ip (k) = 
if i (k) + ip 2 {k) and can therefore conclude from lemma 12.61 that ip' (k) > 0. Plugging 
this in (I3.28P together with x > shows that x' (k) > 0. We thus get that all nodal 
points on the first lead move towards the cycle, as k increases. A similar derivation 
leads to the same conclusion for the nodal points on the second lead. The second real 
contra-phasal solution is a reflection of the one mentioned above and therefore its nodal 
points obviously also move towards the cycle. 

□ 

Lemma 3.10. Let k be a value at which a nodal point is positioned on a vertex ofT. 
The following scenarios exist for the dynamics of the mentioned nodal point. 

(1) The nodal point had arrived to the vertex from a lead. Then, upon entering the 
cycle the nodal point will either split into two nodal points or merge with another 
nodal point arriving from the cycle. The set of k values at which these events 
happen is {k*} := { 2 b+ 2c P} ■ The split events happen at k e {kl, k%, k^, ...} 
and the merge events at k G {k 2 , k\, k$, ...}. 

(2) No nodal point arrives to the vertex from the lead during this event. The nodal 
point had therefore arrived to the vertex from the cycle. It will just flow to the 
other edge of the cycle. These events happen at k values for which 7 (k) G 7rZ. 

Proof. When a nodal point enters the cycle from one of the leads, say the first one, 
fi(k*; 0) = 0, and we have from property 13.71 that on the second lead -£^f2(k*; 0) = 0. 
We therefore have that the restriction of / to the cycle during such an event is equal to 
an eigenfunction of a single edge of length 2b + 2c with Dirichlet vertex conditions at 
its endpoints. This implies that the entrance events occur at k* = 2 b+ 2c P- These events 
are of two types (explanation follows): 

(1) At A?2 m _ 1 = 2b l 2c (2m — 1) the entering nodal point splits into two new nodal 
points which continue to move in the cycle. Hence the total number of nodal 
points increases by one. 

(2) At k 2m = 2 b +2 2m the entering nodal point merges with another nodal point 
coming towards it from the cycle. Hence the total number of nodal points 
decreases by one (see figure 13. 5p . 
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During an entrance event, k = k*, the nodal point is positioned on a vertex of F and 
about to enter the cycle. We observe that the number of nodal points on the cycle 
must be even. This implies that at the entrance event the mentioned nodal point either 
merges with another nodal point form the cycle (so that the number of nodal points 
on the cycle remains unchanged), or splits into two nodal points (which increases this 
number by two). The occurrence of a split or a merge event is determined by the 
values of f(k*; •) restricted on the cycle. As mentioned before, this restriction is an 
eigenfunction on the edge of length 2b + 2c and it is therefore equals sin ( 2 b+2c x ) U P ^° 
a multiplicative scalar. For an even value of p, this function has opposite signs in the 
vicinity of the endpoints of the edge. This means that when the nodal point is located 
exactly on the vertex of T, the two nodal domains of f(k*\ •) on the cycle which are 
bounded by this nodal point have opposite signs, (see figure 13.51 -during) . 



2b H 



before 





after 




FIGURE 3.5. A description of a merge event. The nodal points (2a) 
and (2b) merge and become the nodal point (2). The signs of the nodal 
domains of f^ k ' are marked with squares. 

However, a short while before this event, the neighborhood of this vertex was con- 
tained in a single nodal domain with a definite sign. The /^-continuity of the solution 
implies that this is possible only if a short while before the event there was another 
nodal point in the vicinity of the vertex that has disappeared while merging with the 
nodal point at the vertex (see figure 13.51 -before) . A similar reasoning shows that split 
events occur for odd p values. 

We have treated by now the possibility that the nodal point at the vertex had arrived 
from the lead. It might also happen that f(k; •) equals zero at a vertex of T when 
f(k; •) vanishes on the lead which is connected to that vertex. For the real contra- 
phasal solution given in (I3.26P this happens exactly at 7 (k) = 7rZ. This event would 
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happen only on vertex number two for that solution (and on vertex number one for 
the reflected solution). These events do not change the number of nodal points on the 
graph, and therefore we do not need to keep track of them. 

□ 

Lemma 3.11. Let a, b, c be positive real numbers such that - c ^ Q and Q and T, T 
be the graphs described above. The number of nodal points on T of a real contra-phasal 
solution on T is increased by one at k such that k 2 G o (r) . 

Proof. When k equals an eigenvalue of T, the solution restricted on T equals an eigen- 
function of T, i.e., either f\(k; a) = or f 2 (k; a) = 0. A nodal point is therefore 
positioned on the boundary of T, and from lemma [3791 we deduce that this nodal point 
moves towards the cycle, increasing by one the number of nodal points on T. It is only 
left to verify that there is no simultaneous split or merge events which further change 
the total number of nodal points. Namely, we show that {k*} and o~ (T) are disjoint 
sets. Assume the contrary: k* G a (r) for some p. By definition, fj(k*; 0) = for 
j G {1, 2}. Assume without loss of generality that fi{k*; 0) = 



or a) = 0. If f^k, 



Then, since k* G a (T) 
then k* G \7L 



we also have that either /i(/c*; a) 

and applying lemma [3.101 gives 2fe T 2e G ^Z, which contradicts the incommensurability 
assumption. Otherwise, if f[(k*; a) = 0, we similarly obtain 26 T 2c e la + a^> an< ^ a g arn 
get a contradiction. 

□ 

Lemma 3.12. Let the set {k^} , as defined in lemma X3.1(\ be the set ofk values at which 
merge and split events occur, and let k^ = 0. Denote d p := \\k 2 E a (T) | k*_ 1 < k < k* } | 
, the number of eigenvalues of T that occurred between two consequent merge/split 
events. Then 



(3.29) 



b + c 



V 



b + c 



(p 



+ 1. 



Proof. The following two observations concern the set D\ (k) U Di (k), which gives the 
positions of the nodal points on the leads. 
The spectrum of T may be characterized as 

(3.30) a(T) = {k 2 \ae D 1 (k)UD 2 (k)} . 
The merge/split events happen at 

(3.31) {k* p } = {k \ 0eD 1 (k)UD 2 (k)}. 

We denote Q (k) := (0, a] D {Di (k) U D 2 (k)} and describe how it changes with k. 
Lemma [3.91 implies that the values of Q (k) continuously decrease with k. In addition, 
the first observation gives that \Q (k) \ increases at k G cr (T), when a nodal point enters 
r. The second observation shows that \Q (k)\ decreases by one at k G {kp} , when a 
nodal point enters the cycle. It is therefore evident that during the interval (k^_ 1 , /c*] , 
\Q (k) \ decreased a single time (at k*), and the number of times it increased is given by 
d p , the number of eigenvalues in this interval. We conclude that 

(3.32) d p =\Q{k;)\-\Q(k* p _ l )\+l. 

It is easy to see that Di (k*) U D 2 (k*) = ^Z, and therefore 

\Q(K)\ 







2k* p a 


(0,al n -?-z 
v ' J 2k; 








TT 
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Substituting k* 



26+2c 



p (lemma T3.1 Op and plugging this in (I3.32j) gives (13. 29 p . 



□ 



We now have all the required information to obtain an expression for /i n , the number 
of nodal points on T at k = k n . 



Proof. [Proof of theorem 13.3 



In order to prove (13 . 18[) we need to keep track of all the events which affect the number 
of nodal points on the graph. These include the eigenvalues of the original graph, {k n }, 
and the merge/split events, {k*}- Fi gure [3761 shows a possible scenario for such a stream 
of events. In this figure, the value of /i n — n is shown for each eigenvalue. The bounds 
on n n — n can be obtained from (11. 1 U j) with a slight modification due to the additional 
nodal point positioned on the Dirichlet vertex of the graph: 0< fi n — n < 1. The value 
of n n — n differs from /i n _i — (n — 1) if and only if a merge/split event occurred in 
between the corresponding eigenvalues. We therefore conclude that the value of fi n — n 
depends on the parity of the number of merge/split events that occurred before k n . 



k* 

I 1 h — (§) h- 



k* k* 
2 k 4 k 5 k 6 3 

ImI- 



li n -n 

di di = 2 



1 

da = l 



H 1 — h- — (S) h 



4 = 3 



Figure 3.6. An example of a stream of events needed to keep track of 
the number of nodal points. 



Namely, 



/i n -n — mod 2 p, 

where p is an integer such that 

kp < k n < k p ^. 
By the definition of d p (see lemma I3.12p this is equivalent to 

p p+i 
< n <^di, 

i=\ i=l 



which by (I3.29P evaluates to 
a 



-P 



+ p < n < 



_b + c 

Since n, p are integers and ^ ^ Q, 



b + c 



ip + K 



Multiplying through by a ^^ c we get 



+ (p + 1) • 



T^—p + p<n< — ^— (p + 1) + (p + 1) . 
+ C + c 



b + c , _ 

p < ; n < (p+1) 

a+b+c K ' 
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and conclude that 



P 



b + c 



-11 



-n 



a + b + c 

The number of nodal points on the graph is therefore given by 

b + c 
a + b + c 

We now wish to turn this into a formula for the nodal count, u n . The relation between 
fi n and v n depends on whether the n th eigenfunction has nodal points on the loop as 
demonstrated in figure 13.71 



fjL n = n + mod 2 



A*n = 4 u n = 3 

(b) 



FIGURE 3.7. The two possible relations between the number of nodal 
points, fi n , and the number of nodal domains, v n . 

If it does have nodal points on the loop then u n = fi n — 1 (figure 1577T b)). and in the 
case it does not, v n = fi n (figure l3T7T a)). 

We therefore have that for the first d\ = \_-^\ + 1 eigenvalues (when there are still no 
nodal points on the loop) the nodal count is 

b + c 



v n = n + mod 2 



-n 



n. 



a + b + c 

where the second equality is due to n < + 1. For the rest of the nodal count, 

b + c 



n — 1 + mod 2 



a + b + c 



-n 



□ 



4. The nodal count of graphs with disjoint cycles 

4.1. Graphs with — 1: a dynamical approach. In this section we will discuss 
the nodal dynamics on a graph with one cycle (i.e. /3 — 1) and a lead attached to a 
general position on the cycle, see figure 14.11 The discussion will not be formal, as we 
will prove the results by other methods in section 14.21 

We have seen in section I3TT1 that, as k increases the nodal points (zeroes) travel along 
the lead in the direction of the graph. Consider the quantity 5 n = fi n — n + 1. This is 
the "surplus" of zeros due to the graph not being a tree. Bound fll.lOp implies that 5 n 
can be equal to either or 1. The change in this quantity from eigenvalue k = k n -\ to 
eigenvalue k = k n can be attributed to the following three causes: 

(1) The increase in the index n (the change in d is — 1). 

(2) A zero entering the graph from the lead. Upon entering, the zero can either 
merge (M) with a zero already present on the cycle or split (S) into two zeros. 

(3) A zero entering a tree. This zero can either split off a zero traveling on the cycle 
or it can be a result of two zeros from the cycle merging together. 



NODAL COUNT ON A FAMILY OF QUANTUM GRAPHS 



25 




[ lead 



Figure 4.1. An example of a graph T with (3 = 1. We view it as a cycle 
with trees grafted upon it. 

Another notable event is a zero passing through a vertex where a tree is attached. 
We did not list it above since an event of this type does not affect the nodal count. 
Similarly, when a zero is traveling through the tree, we know (see [T5| [T7J [2D]) that the 
number of zeros does not change. 

As already explained in section I3TT] event (2) happens exactly once between each pair 
of eigenvalues fc n _i and k n , since an eigenvalue corresponds to the Neumann condition 
and the entrance event corresponds to the Dirichlet condition satisfied at the attachment 
point. If event (2) is a split, the contribution to 5 is +2, otherwise it is 0. However, if 
we consider the total contribution of events (1) and (2), we get +1 from a split and — 1 
from a merge. This is the same as a contribution of a type (3) event, when the split 
results in +1 (number of zeros on the cycle stays the same but another zero appears on 
a tree) and the merge in —1 (the number of zeros on the cycle reduces by 2, while one 
zero enters a tree). 

The first eigenf unction has constant sign, so 5i = and no events happen until 
k = k\. Since the contribution of type (1) is now absorbed in the contributions of type 
(2), the value of 5 n is the total number of splits minus the total number of merges up 
to k = k n . On the other hand, 5 n is restricted by the nodal bound to be either 1 or 0, 
therefore it is equal to the parity of the total number of S/M events. 

There are exactly n — 1 events of type (2) happening until k = k n . To count the 
number of events of type (3), we consider an auxiliary graph V®, obtained from T by 
removing all edges belonging to the cycle and imposing Dirichlet conditions on the 
points where the trees were connected to the cycle. The graph is a collection of trees 
that were grafted on the cycle. Since a zero entering a tree signals that the Dirichlet 
condition is satisfied on the tree, the corresponding value of k is in fact an eigenvalu^l 
of r<ft. And the number of events of type (3) is thus equal to the number of eigenvalues 
of Tg that are smaller than k n . To summarize, 

S n = mod 2 (n - 1 + N®(k n )), 

where N%(k n ) is the spectral counting function of the graph T$. Thus we can fully 
predict the nodal count using the spectra of two graphs, V and IV The discussion 
above captures the dynamics of the zeros, but it is relatively difficult to formalize. 
Instead we will prove the formula for /x n by other methods, which, although not very 
pictorial, allow us to extend the argument to the case of non-zero potential V(x). 



The corresponding eigenfunction is identically zero on all trees apart from the one with Dirichlet 
condition satisfied. 
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4.2. Graphs with /3 = 1: a formal proof. In this section we prove the formula that 
was informally derived in section 14.11 

Theorem 4.1. Consider the Schrddinger operator (II. 2p on a connected graph with a 
single cycle. Let the n-th eigenvalue be simple and the corresponding eigenfunction be 
non-zero on the vertices. Then 

(4.1) /i n = n-l + mod 2 (n-l + A^ (A„)), 

where N$ is the spectral counting function of the disconnected graph obtained by remov- 
ing the cycle and putting Dirichlet conditions on the new vertices. 

Proof. From the nodal bound for graphs with one cycle (i.e. with (3 = 1) we know that 
fi n is equal to n — 1 or n. The first step of the proof is to observe that the number of 
zeros on the edges that do not belong to the cycle is equal to N$(\ n ). We will prove 
this statement below. Then n n — N$(X n ) is the number of zeros on the cycle, and has 
to be even, as explained in section 11.21 

First, assume that « n = n. Then, the quantity 

H n - 1 - iV (A n ) =n-l- iV (A n ) 

is odd and therefore 

mod 2 (n - 1 + AT (A„)) = mod 2 (n - 1 - iV (A n )) = 1, 
where we used the fact that 

mod 2 (a + b) = mod 2 (a — b) 

for any integer a and b. Thus the right-hand side of equation (14. ip evaluates to n which 
is the right answer. 

If fi n is not equal to n, then it is equal to n — 1 and we have 

mod 2 (n - 1 + iV (A n )) = mod 2 (ji n - N^(\ n )) = 0, 

since u„ — A^(A n ) is even. Thus equation (14. ip still holds. 

Now we prove that iV0(A„) is indeed the number of zeros on the subtrees of the 
graph. To shorten the formulas we introduce the following notation. We denote the 
n-th eigenvalue by A and the corresponding eigenfunction by F. Break up the original 
graph into the cycle and the trees Tj. For each tree Tj we choose as a root the vertex 
that was its contact point with the cycle. We can ensure that each root has degree 1: 
if necessary we can split trees that share a root. On each tree the vertex conditions are 
inherited from the graph, but we still need to specify the conditions on the root. We 
will consider two versions of each tree. The first, Tj t p has the condition on the root r 
chosen to be satisfied by the function F, restricted to the tree. That is we chose the 
constant a in the 5-type condition to be a = F'(r)/F(r). The second version of the 
tree, denoted Tj j00 , has the Dirichlet condition on the root. 

Denote by Too, the disjoint union of the graphs Tj jOC . We observe that 

j 

Thus we only need to prove that iVy. x (A) gives the number of zeros of F on the subtree 
Tj j00 . Since, by construction, the restriction of F is an eigenfunction of Tj t p with the 
eigenvalue A, we have A = \ m (Tj^p) for some m. By the strict interlacing, theorem lA.il 

A m -i(Tj >00 ) < \ m (Tj t p) < A m (T, j00 ), 



NODAL COUNT ON A FAMILY OF QUANTUM GRAPHS 



27 



and, therefore, iV^ (A) = m — 1. On the other hand, the nodal count on trees, 
equation fll.lOp . gives u m (T., )0O ) = m — 1. This concludes the proof. □ 

4.3. Number of zeros on a graph with disjoint cycles. In fact, the formula of 
the previous section can be extended to > 1 as long as the cycles do not share any 
vertices. 

Theorem 4.2. For a connected graph containing j3 disjoint cycles, let the n-th eigen- 
value be simple and the corresponding eigenfunction be non-zero on the vertices. Then 

H n = n - 1 + mod 2 (n - 1 + A^A^)) , 

3=1 

where Nj$ is the spectral counting function of the disconnected graph obtained by re- 
moving the j-th cycle and putting Dirichlet conditions on the new vertices. 

Proof. Denote the n-th eigenvalue by A and the corresponding eigenfunction by F . 
Choose an arbitrary cycle and let e\, . . . , be the edges incident to it. Since the 
cycles are disjoint, these edges do not belong to any cycle. Choose points C\, . . . , c&, 
one on each edge, so that the function F is non-zero at these points. If the graph is cut 
at these points, we obtain k + 1 disjoint subgraphs, Tj, j = 0, . . . , k (the 0-th subgraph 
contains the chosen cycle), see figure |4~21 




Figure 4.2. An example of a graph with 4 disjoint cycles. In part (a), 
the middle cycle is chosen and the cut-points are labeled. In part (b) the 
graph is split up into subgraphs. The central part To contains exactly 
one cycle, while the other parts can contain more or fewer. 



Define 



F(cj 



-3) 

where the derivative is taken away from the chosen loop. We impose 5-type conditions 
on the newly formed vertices. The vertex Cj belonging to Tj will get the condition with 
coefficient a c . and its counterpart belonging to the subgraph Tq will get the condition 
with coefficient —ce c .. This way, the appropriately cut function F is still an eigenfunc- 
tion on all k + 1 subgraphs and A is the corresponding eigenvalue. This allows us to 
define rrij by 

A = A m ,.(r i ). 
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Lemma 4.3. The numbers rrij are well- defined and satisfy 

k 

(4.2) 71-1 = ^(77^-1). 

3=0 

Proof of the lemma. Let r c denote the disjoint union of the graphs Tj, j = 0, . . . , k. 
First of all, we apply theorem IA.2I k times (for k cuts) to obtain inequalities 

A n -i(r c ) < V-i(r) and A n+ i(r) < A n+fe+ i(r c ) 

On the other hand, simplicity of the eigenvalue A n (r) = A means that A„_i(r) < A < 
A n+ i(r) and therefore 

(4.3) A n _x(r c ) < A < A n+fc+1 (r c ). 

Finally, out of F we can form at least k + 1 linearly independent eigenfunctions of the 
graph r c : functions that are restrictions of F on one of the parts Tj and identically 
zero on all the others. All these eigenfunctions have eigenvalue A. Combining this 
observation with inequality (14. 3 P we conclude that A has degeneracy exactly k + 1 in 
the spectrum of r c and therefore is a simple eigenvalue of every part Tj. Thus the 
numbers uij are well-defined. 

Finally, since the spectrum of T c is the superposition of spectra of Tj , equation (14. 2 h 
is equivalent to the statement "there are n — 1 eigenvalues of T c that are strictly less 
than A", which is also obvious from inequality (14. 3 P and the fact that A n (r c ) = A. □ 



(a) (b) 

V i?i \ I R 




■2 

D D 



D D 

R4 ' s Rs 



FIGURE 4.3. The subgraph r and the graph R obtained after removing 
the cycle (it appears shaded on part (b) of the figure). The graph R is a 
disjoint union of four subgraphs Rj. 

We now want to use theorem 14. II to find the number of zeros of the function F on the 
graph Tq. Let R be the graph obtained from Tq by removing the cycle and imposing 
Dirichlet conditions on the new vertices. This graph is a disjoint union of the graphs 
Rj, j = 1, . . . , k, see Fig. 14.31 Therefore, we have 

k 

(4.4) N r (A) = ^2Nr.(A). 

3=1 

According to theorem 14.11 the number of zeros of F on the subgraph T is 

/i(A, r ) = m - 1 + mod 2 (m - 1 + N R (A)) . 
Extracting m — 1 from equation (14. 2 p and using equation (14. 4p we get 

(4.5) /i(A, r ) = m - 1 + mod 2 \ n - 1 + ^ (rrij - 1 + N Rj (A)) 
where we used 

mod 2 (a + b) = mod 2 (a — b) 
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for integer a and b to change some signs. Define now the graph G by removing the 
chosen cycle from the original graph T and imposing the Dirichlet conditions on the 
new vertices. Similarly to the graph R, the graph G is a disjoint union of k subgraphs 
Gj, see Fig. 14.41 and 



'31 

k 

(4.6) N G (A) = J2N Gj (A). 

i=i 

If we were to cut the graph Gj at the point Cj the two parts would be exactly Tj and 
Rj. This suggests the following lemma. 




FIGURE 4.4. The original graph T and the graph G obtained by removing 
the chosen cycle. The graph G is a disjoint union of graphs Gj. 



Lemma 4.4. For every j = 1, . . . , k, 

(4.7) N Gj {A) = mj -l + N Rj {A). 

Proof of the lemma. First we observe that A belongs to the spectrum of the graph 
Tj and does not belong to the spectrum of Rj or Gj. This can be shown by the strict 
interlacing, theorem lA.il applied to the graph T (corresp. Go) by changing the condition 
from Neumann to Dirichlet at the vertex where Gj (corresp. Rj) was connected to the 
cycle. 

Denote by Gj jC the disjoint union of the graphs Tj and Rj. Let integer q be such that 

A = A,(G iiC ). 

Since A = A m .(r.,), we have that q = rrij + Nr. (A). On the other hand, by theorem 

\ q „ l (Gj)<A = \ q (G jiC )<\ q (G j ) 
Therefore, N Gj (A) = q — 1 which concludes the proof. □ 

Combining equations (14. 6 j) and (14. 7p with formula ( 14. 5 P we arrive to 

(4.8) //(A, r ) = m - 1 + mod 2 (n - 1 + N G (A)) . 

We should also note that N G = Nj t $, where j is the number of the cycle that was chosen. 

Since we chose an arbitrary cycle, equation (14.81) is valid for every cycle. The conclu- 
sion of the theorem is just the sum of equations (14. 8 p over all cycles with an application 
of the analogue of equation (14. 2p . □ 
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5. A DISCUSSION 

5.1. An approximate derivation of an exact nodal count formula. We will now 

present an alternative way to get the nodal points count formula f 1 3 . 1 S j) of the graph 
given in figure iSTTT a). The derivation is most appealing, but involves an approximation 
that cannot be justified. We present it here because it makes use of an idea which 
has been used in other contexts. Also we find that an unjustifiable approximation that 
reproduces the exact nodal count formula carries information about the graph in its 
own right. 

(a) (b) 



2b 



N 




D 



FIGURE 5.1. (a) The quantum graph whose nodal count we compute, 
(b) The same graph with its cycle removed (appears shaded). 



We start by rewriting the formula (14. 1|) with a slight modification, due to the nodal 
point which is positioned on the Dirichlet boundary vertex of the graph: 

H n = n + mod 2 (n + N^(k n )) . 

The spectra of the two edges which appear in figure l5TTT b) are {^ n } neN and (n + |) } 
Their spectral counting function is therefore 



net 



ka 




ka 


1 




+ 






TT 




TT 


+ 2. 



N 9 {k) 

Plugging it above and using the identities mod2 ( L X J + |_ x + |J ) = m °d2 ( |_2rcJ ) and 
mod2 (x + y) = mod 2 (x — y) we obtain 

2 k'ri Ob 



fi n = n + mod2 n — 



71 



We can get an approximate expression for k n from the Weyl term of the spectral count- 
ing function of the graph, 



N(k)* 2(a + '' + C h, 



TT 



by its inversion, i.e., 



TT 



" 2(a + b + c) 



n. 



Using this approximation gives 

fx n pa n + mod 2 ( n 

= n + mod2 
which is the exact result, (j3.18p . 



a + b + c 
b + c 



n 



-n 



a + b + c 
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One should note that the last step of the derivation, which involves an approximation 
of k n by inverting the Weyl term, cannot be justified. Moreover, the floor function 
is a discontinuous function and it is therefore expected that an approximation of its 
argument would lead to a completely wrong result for some portion of the sequence. 

From the exactness of the final result, we conclude the following property of the 
spectrum 

2a 7i 
Ti 2(a + b + c)' 

Numeric examination reveals that the equality hold for the arguments of the mod2 as 
well, namely 

(5.1) — k n = — — — ^— — -n 



mo do 





2a 




( 


kn 


J = mod 2 ( 




7T 





-n 



t vr 2(a + b + c) 

The above relation connects the spectrum and the lengths of the graph's edges. Having 
such a relation for our graph makes N%(k n ) expressible in terms of the parameters 
a, b, c and enables to turn the nodal count formula, (14. ip . into a formula which contains 
geometric properties of the graph, rather than spectral ones. In short, the special 
nodal count formula is a direct consequence of a purely spectral identity - a connection 
between the graph's spectrum and the spectral counting function of its subgraphs. 

The novelty of this result makes one wonder to what extent it can be generalized 
to other graphs. Even if such an exact result is not reproduced, one may still use 
approximations of the type above and try to estimate the errors caused by them. 



5.2. Periodic orbits expansions. Wishing to express the nodal count formula A3. 18[) 
as a periodic orbits expansion, we notice that mod2 ( [^J ) is an odd periodic function 
(of period 2), whose Fourier transform is: 



1 °° 2 
mod 2 ( W ) = --y:- rFT ^sm((2* : + l 

k=0 v ; 



71X ) . 



Denoting a := J^_ c , the normalized length of the loop, we can rewrite (13 . as 



following: 



1 oo 2 

/*» = n + o - 22 (2k + 1)tt Sin ^ 2k + X) 7ran) ■ 

fe=0 ^ ' 



We therefore get that the nodal points sequence is expressed in terms of lengths of 
periodic orbits on the graph. One should note that the only periodic orbits that appear 
are odd repetitions of the graph's cycle. They appear with harmonically decaying 
amplitudes. This causes to seek for a more direct derivation of the periodic orbits 
expansion which may also explain the meaning of the amplitudes and the absence of 
other periodic orbits. Furthermore, the formula (14.11) . which holds for any graph with 
a single cycle, may also be turned into an expansion of a similar type. We recall that 
for quantum graphs there exist an exact periodic expansion for the spectral counting 
function. Therefore, the spectral counting function of the subgraph, N$(k n ), can be 
expanded and plugged in formula (14. ip . This would yield an expansion which still 
involves the spectral information, {k n }. Having an approximate inversion of the spectral 
counting function of the whole graph then enables to further get a periodic orbits 
formula which involves only geometric properties of the graph. Such spectral inversion 
attempts were recently carried out with a high degree of success [10, 42J. It is therefore 
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evident that the obtained result leads to a wide field of further questions and open 
research possibilities. 
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Appendix A. Interlacing theorems for quantum graphs 

Eigenvalue interlacing (or bracketing) is a powerful tool in spectral theory. In par- 
ticular, in the graph setting, it allows to estimate eigenvalue of a given graph via the 
eigenvalues of its subgraphs, which may be easier to calculate. Here we quote the 
theorems that are used in the proofs of the formulas of the present manuscript. The 
theorems are quoted in the form they appear in |35j . 

The first theorem deals with choosing a vertex on the graph T and changing the 
parameter a v of the extended 5-type condition at v (see equation fll.4p ). We remind 
the reader that a v = oo corresponds to the Dirichlet condition at the vertex which 
essentially disconnects the edges meeting at the vertex. 

Theorem A.l (Interlacing when changing a parameter). Let T a r be the graph obtained 
from the graph T a by changing the coefficient of the condition at vertex v from a to a' . 
If — oo < a < a' < oo, then 



If the n-th eigenvalue of T a i is simple and the corresponding eigenfunction is nonzero 
on the vertices, the inequalities are strict. 

The second theorem deals with the situation when the graph I~" is obtained from V by 
gluing two vertices together, or, equivalent ly, V is obtained by cutting the graph V at 
a vertex or at a point on an edgejj When gluing the vertices together, their respective 
parameters a VQ and a vi get added. 

Theorem A. 2 (Interlacing when gluing the vertices). Let T be a compact (not necessar- 
ily connected) graph. Let vq and v\ be vertices of the graph V endowed with the 5 -type 
conditions with the parameters «o and ot\ (see definition \1.1\) . Arbitrary self-adjoint 
conditions are allowed at all other vertices ofT. 

Let V be the graph obtained from V by gluing the vertices Vq and v± together into a 
single vertex v, so that £ v = £ VQ U £ V1 , and endowed with the 5-type condition with the 
parameter a v — ao + a%. 

Then the eigenvalues of the two graphs satisfy the inequalities 



In addition, if A n (r') is simple and the corresponding eigenfunction is nonzero on 
vertices and not an eigenfunction of T, the inequalities are strict. 

2 Any point on an edge can be viewed as a vertex of degree 2 
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(A.2) 



A n (r)<A n (r')<A„ +1 (r). 
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An intuitive explanation for the above result is that by gluing vertices we impose an 
additional restriction: the continuity condition. This additional restriction pushes the 
spectrum up. 
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